Molecular Hydrogen Fluorescence in the Eridanus Superbubble 
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ABSTRACT 

The first FUV (AA 1350-1750A) spectral imaging observations of the Eridanus su- 
perbubble, obtained with the SPEAR/FIMS mission, have revealed distinct fluorescent 
emission from molecular hydrogen. In this study, the observed emission features were 
compared with those from a photo-dissociation region model with assumed illuminat- 
ing stellar fields. The result showed rather high line ratios of Ii58o/Ii6io 5 which may 
imply the existence of high-temperature molecular clouds in the region. The H2 fluo- 
rescence intensity showed a proportional correlation with H-a emission, indicating that 
the fluorescence and the recombination emission have similar physical origins. 

Subject headings: ISM: molecules — ISM: supernova remnants — ultraviolet: ISM 

1. Introduction 

The Eridanus region contains a large nearby interstellar superbubble. This region was dis- 
covered as an extended 15° X-ray 'hot spot' (Williamson et al. 1974; Naranan et al. 1976). Sub- 
sequently, the Eridanus Loop was identified as an expanding shell of H I by Heiles (1976) and 
was extensively imaged in H-a emission (Reynolds et al. 1998; Boumis et al. 2001), revealing a 
large filamentary shell structure. Reynolds and Ogden (1979) argued that Barnard's Loop and 
the Eridanus shell were both created by one or more supernovae in the Orion OB association and 
were maintained by the intense stellar winds and ionizing radiation from those stars. ROSAT ob- 
servations (Snowden et al. 1995) support the idea that the huge interstellar cavity is filled with 
several-million-degree gas. Surveys of CO (Magnani, Blitz, and Mundy 1985; Aoyama et al. 2002) 
show that the cavity is bounded in part by molecular gas. The CO clouds generally coincide with 
H I filaments; however, on a smaller scale, the clouds are located where H I is weak within the 
filament. In the FUV spectral range, early (Murthy et al. 1993) and improved (Kregenow et al. 
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2006) spectral observations of diffuse emission from the Eridanus Loop over limited angular fields 
show that gas at intermediate temperatures, T~ 10 5 K, exists between the hot superbubble interior 
and the surrounding material. 

We present the first FUV spectral imaging observations of the entire Eridanus Loop. The data 
were obtained with the Spectroscopy of Plasma Emission from Astrophysical Radiation (SPEAR), 
also known as the Far-ultraviolet Imaging Spectrograph (FIMS). These SPEAR data show H2 
fluorescent emission throughout the Eridanus Loop field and provide new insights into the geometry 
and physical conditions of the molecular environment around this superbubble. 

2. Observations and Data Reduction 

The Eridanus region was observed with SPEAR, the main payload of the Korean satellite, 
STSAT-1, which was launched on September 27, 2003. SPEAR is a dual-channel FUV imaging 
spectrograph (S channel 900-1150A, L channel 1350-1 750A, A/AA ~ 550, where AA refers to the 
half energy width) with a large field of view (S: 4.0° x 4.6', L: 7.5° x 4.3', angular resolution~ 10') 
optimized for the observation of diffuse emissions. See Edelstein et al. (2006) for an overview of 
the instrument, mission, and data analysis. 

The region bounded by a = 2 h 30 m and 5 h 30 m , and S = +13° and —23°, covering most of the 
constellations Orion and Eridanus, was observed using more than 16 ks of exposure taken between 
26 December 2003 and 26 January 2004. The observation geometry is shown in Figure 1, along 
with an H-a survey map. The H-a features most prominent around the Eridanus supershell are 
Barnard's Loop, the I Orion O-association molecular cloud complex, the H II region surrounding 
A Ori centered at 5 h 32 m , +9°50 m , and H-a filamentary structures. Barnard's loop and the A Ori 
region were excluded from the SPEAR observation to avoid excessive count rates. For the ensuing 
data analysis, the region was divided, based on the H-a emission features, into subregions labeled 
as 1 through 6 (See Figure 1). 

We identified bright stars in our observations by examining an FUV sky intensity map, created 
by accumulating photons and exposure into 2' sky bins. Stars were identified by matching features, 
within 15', to the position of UV bright stars in the Tycho-2 catalog (Wright et al. 2003). Photons 
within 30' of the identified stars were removed from the data set to obtain the residual diffuse 
emission and to minimize stellar contamination. After excluding the contributions of UV bright 
stars, 3.0 x 10 5 of 8.8 x 10 5 detected photons remained. These photons were accumulated, for each 
subregion, into spectra that were binned to 2A, shown in Figure 2(a)-(f), along with the noise levels 
per spectral bin. An inherent detector background of ~ 0.013 counts/s/A was measured during 
interleaved shutter-closed periods and was subtracted from these spectra. 
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3. Spectrum Analysis 

The Eridanus 1350A-1630A spectra show emission features that correspond to the H2 fluores- 
cence predictions of Sternberg (1989), notably the prominent bands near 1605A. Other prominent 
features were identified with atomic lines such as C IV at 1550A and Si II at 1530A (See Kregenow 
et al. (2006) for a discussion of atomic emission in the Eridanus Loop.). To investigate the physical 
conditions of the H2, a photo dissociation region (PDR) radiation code (Black and van Dishoeck 
1987) was used to generate a model FUV spectrum. The continuum level was taken to be a linear 
function of wavelength, and the sensitivity of the fitting while varying the slope and y-intercept 
of the continuum was investigated based on the x 2 minimization method. The fit models for each 
subregion are shown in Figure 2. 

There are several physical parameters which affect the resultant spectrum: cloud density, 
excitation temperature, total hydrogen molecule column density, and radiation field; we varied 
only the excitation temperature in our model with other parameters fixed. For example, ultraviolet 
fluorescent emission of H2 is produced following the absorption of photons in the Lyman (B x £+ 
- X 4 X+) and Werner (C 1 IT M - X X S+) bands. Hence, the stars in the I Orion OB association 
were considered as the likely and dominant source of H2 emission. The total radiation field from 
the most luminous stars of the I Orion association (Reynolds and Ogden 1979) accounts for the 
presumed radiation field scale factor (Habing 1968) of x ~ 1-9, assuming that the distance from the 
association is ~ 200 pc and there is no extinction. Heiles, Haffner, and Reynolds (1999) estimated 
the electron volume density, n e , of the ionized gas in the Eridanus superbubble to be 0.8 cm~ 3 , 
which corresponds to the thermal pressure, p/k ~ 1.5 x 10 4 cm~ 3 K, in the ionized gas. The 
density of the neutral gas was adjusted to have the equilibrium thermal pressure, p/k ~ 1.5 x 10 4 
cm" 3 K, according to the excitation temperature, which is used as a parameter in our model. The 
fractional abundance of H2, defined as (/) = 2(n(H2))/[2(n(H2)) + (n(H I))], was assumed to have 
an average value of 0.17, as obtained from the Copernicus observations (Savage et al. 1977). The 
color temperature of the incident radiation was assumed to be 3 x 10 4 K. 

With these parameters fixed, the molecular cloud excitation temperature was varied to inves- 
tigate the sensitivity of the model spectrum compared to the observed data. The model spectrum 
was scaled to match the average observed intensity over the 1350-1650A spectral region, where H2 
emission was apparent in both the model and the data. The 1356A airglow line did not appear in 
the observed spectrum. The sensitivity of the x 2 value to the parameters comes from the relative 
intensity variations of each vibrational transition band, which are governed by the distribution of 
H2 among the vibrational excitation states in the electronic ground state (X 1 S+). For example, 
the intensity ratio, Ii58o/Ii6io> increases monotonically from 0.65 to 0.8 as the excitation tempera- 
ture increases to 1,000 K, according to our model. Thus, a unique excitation temperature can be 
determined if the line ratios are measured. 

The results for each subregion, the maximum fluorescent emission intensity, ImaxH 2 > the con- 
tinuum intensity I CO nt, and the intensity ratio, Ii58o/Ii6io> are listed in Table 1. The model spectra 
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with T ex = 50 K and T ex = 1,000 K are plotted for Region 4 in Figure 2(g), with the observed 
spectrum. As can be seen in the figure, most of the lines in the spectral range are identified as 
Lyman bound-bound transitions (BX v' — v"), where v' and v" represent the vibrational excitation 
of a hydrogen molecule at each respective excitation state. The observed spectrum best fits with 
a cloud excitation temperature of 1,000 K, with xl = 2.2, though the signal-to-noise ratio of our 
observation was not sufficient to discriminate the excitation temperature accurately. 

4. Discussion 

Our detection of H2 fluorescence in Eridanus is consistent with previous observations in other 
locations in that fluorescence is found where H2 is exposed to a strong UV radiation field. The 
correlations among the subregions' fluorescent emission intensity, the H-a emission intensity, and 
the color excess value are shown in Figure 3. The H2 fluorescence intensity shows a strong correlation 
with both the H-a emission intensity and the color excess. This suggests that the H-a emission 
and hydrogen molecular fluorescence have similar physical origins and are governed globally by a 
strong and similar UV radiation field, likely arising from the I Orion OB association. We expect 
that the H2 abundance is proportional to that of both hydrogen atoms and dust particles, since 
dust particles are known as the main birth place of H2 in the interstellar medium. 

The most conspicuous features in the H-a emission map shown in Figure 1, are the filamentary 
structures elongated along a = 4 h (Arc A; Region 3) and a = 3 h 20 m (Arc B; Region 5). It appears 
from Figure 3(a) and (b), and from Figure 1, that UV photons from the Orion OB association 
undergo extinction by the Arc A and Arc B filaments, since Region 6 shows relatively low H-a and 
H2 fluorescence, while having a color excess value comparable to those of Region 1, 2, and 5. The 
H2 fluorescent emission from Region 5 (Arc B) is relatively high among regions with comparable 
dust extinction (Regions 2 and 6), and it is similar to that of Region 4, even though the dust 
extinction for Region 4 is higher and the line-of-sight distance to the I Orion OB FUV sources is 
likely smaller. This observation can be explained if the gas in Region 3 (Arc A) is not in the line 
of sight that extends from the radiation source (the I Orion OB association) to the gas in Region 5 
(Arc B), i.e., the distances to Arcs A and B are different and the Eridanus shell structure is open 
from the I Orion to Arc B. This geometry is consistent with the considerations of Boumis et al. 
(2001) and Welsh, Sallmen, and Jelinsky (2005) in that both studies argue that Arcs A and B are 
part of a complex of individual shells viewed along the same line of sight but at different distances 
(Arc A >500pc, Arc B ~150pc). The observational results do not rule out Reynolds and Ogden's 
(1979) description, in that the both Arc A and Arc B can still be exposed to the radiation from 
the I Orion OB association that causes the H2 fluorescent emission. 

The FUV continuum intensity does not show a strong correlation with H-a emission or H2 
fluorescent emission intensities (See Table 1). Instead, the continuum intensity seems closely tied to 
the geometry of the region. The strong continuum of Region 1 can be attributed to the proximity of 
Region 1 to the radiation source. Region 3 is remarkable in that it shows both a strong continuum 
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intensity and a strong H2 fluorescent emission intensity. 

The observationally derived shell expansion velocity of the Eridanus region, of between 15 and 
23 km s" 1 (Heiles 1976; Reynolds and Ogden 1979) is consistent with the survival of the H 2 in the 
Eridanus region. The observed I1580/I16IO value suggests the existence of high-temperature molec- 
ular gas in the region. According to Draine, Roberge, and Dalgarno (1983), when an interstellar 
shock with v s = 25 km/s is running into molecular gas with a density of n = 1 x 10 2 cm -2 , the 
region in which the gas temperature exceeds 200 K will extend only 2 x 10 16 cm, implying that 
interstellar shocks are not good at producing large column densities of warm gas. An alternative 
explanation, for the heating of molecular gas, might be selective UV excitation by UV lines radiated 
by the hot gas in the bubble, or it may be that other heating mechanisms are able to raise the gas 
temperature where UV pumping of the H2 is taking place. Further studies are needed to clarify the 
contribution and role of UV radiation, supernova shocks, and other possible heating mechanisms. 

In summary, we have presented the first detection of H2 fluorescence in the Eridanus region. 
The fluorescence emission features, with a PDR spectrum prediction model, imply the possible 
existence of molecular gas with an excitation temperature as high as 1,000 K, which is about an 
order of magnitude higher than excitation temperatures generally found for molecular gases in the 
Galactic disk (Savage et al. 1977). We note that high-excitation temperature H2 has also been 
detected in high-velocity shocked gas associated with the Monoceros Loop (Welsh, Rachford, and 
Tumlinson 2002) and in a galactic translucent cloud (Snow et al. 2000). The observed spectrum 
of the H2 fluorescence in the Eridanus region is different in that a substantial column density of 
H2 (~ 1.0 x 10 20 cm" 2 ) is required to produce the observed fluorescence. Thus, most of the H 2 
undergoing photon pumping must have a high rotational temperature, if collisional excitation is to 
account for the observed fluorescent feature ratios. A more rigorous observation and model study 
are expected to reveal the physical conditions of the molecular gas in the Eridanus supershell. 

SPEAR/FIMS is a joint project of KASSI, KAIST (Korea) and U.C., Berkeley (USA) and is 
funded by the Korea MOST and NASA (Grant no. NAG5-5355). 
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Fig. 1. — Observation geometry of the Eridanus region marked over an H-a map. The marked 
rectangular region was observed with SPEAR. The observed region was divided into subregions, 
labeled 1 through 6, for spectral analysis and interpretation. The thick solid line, in the upper 
right corner, designates the SPEAR field of view. 
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Fig. 2. — SPEAR FUV diffuse spectra of Eridanus subregions 1-6, shown in panels (a) - (f). The 
histogram is the observed spectrum. The thick solid line and the dashed line are the best-fit model 
of H2 fluorescent emission and continuum intensities, respectively. The noise levels are shown as 
vertical bars, (g) The continuum subtracted diffuse spectrum of Region 4 with different H2 model 
parameters. The dashed line represents for T ex = 50 K. The solid line represents for T ex = 1,000 
K. The legends, marked as BX v' — v" , show the positions of Lyman bound-bound transitions. 
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Fig. 3. — Comparison of the H2 fluorescent emission intensity (asterisks) for the subregions with 
(a) the H-a emission intensity and (b) the color excesses (Finkbeiner 2002). The measurement 
error ranges are represented with error bars. The solid lines are the best linear fits, and the dashed 
lines show the 1 a uncertainty in b, in the linear fit y = ax + b. 
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Table 1: Observational results for regions in Eridanus. 



Region 


T a 

L cont. 


Imax.H2 




Il58o/Il610° 


H-a(R) d 


E(B-V) c 


1 


1590 ± 165 


825 ± 130 


2.1 


0.8 ±0.4 


5.7 ±0.8 


0.11 


2 


1065 ± 110 


625 ± 90 


2.3 


0.5 ±0.3 


7.0 ±0.9 


0.08 


3 


1555 ± 160 


1430 ± 175 


2.0 


1.2 ±0.35 


17.7 ±2.0 


0.19 


4 


915 ± 100 


1205 ± 155 


2.2 


0.9 ±0.25 


8.0 ± 1.0 


0.14 


5 


720 ± 80 


1190 ± 145 


3.8 


0.85 ±0.25 


9.0 ± 1.2 


0.08 


6 


520 ± 60 


470 ± 80 


2.9 


1.2 ±0.6 


1.3 ±0.4 


0.08 



"'^The intensity of the continuum at 1525A, l CO nt, and the maximum H2 fluorescent emission from 1600 to 1620A, 
Imax.H 2 , in units of photons/cm 2 /sr/A/s. 

c The line ratio, Iisso/Ii6io, is estimated by summing 3 spectral bins (6 A interval), centered at 1580A and 1610A. 
d The H-q intensity values from Reynolds et al. (1998) 
e The color excess values, E(B-V) from Finkbeiner (2002). 



